Robustness of edge states in graphene quantum dots 



M. Wimmer and A. R. Akhmerov 
Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 

F. Guinea 

Instituto de Ciencia de Materiales de Madrid, CSIC, 
Sor Juana Ines de la Cruz 3, E28049 Madrid, Spain 
(Dated: March, 2010) 

We analyze the single particle states at the edges of disordered graphene quantum dots. We show 
that generic graphene quantum dots support a number of edge states proportional to circumference 
of the dot over the lattice constant. Our analytical theory agrees well with numerical simulations. 
Perturbations breaking electron-hole symmetry like next-nearest neighbor hopping or edge impuri- 
ties shift the edge states away from zero energy but do not change their total amount. We discuss 
the possibility of detecting the edge states in an antidot array and provide an upper bound on the 
magnetic moment of a graphene dot. 



I. INTRODUCTION 

The experimental discover};^ ^ of graphene, a mono- 
layer of carbon atoms, has opened room for new elec- 
tronic devices (for reviews, see Refs. 3-5). A peculiar- 
ity of finite graphene sheets is the existence of electronic 
states localized at the boundary, so-called edge states. 

A crystallographically clean zigzag edge was theoret- 
ically predicted to sustain zero-energy edge state^^^^. 
Later, it was shown^ that any generic graphene boundary 
not breaking electron-hole subband (sublattice) symme- 
try also supports these zero energy edge states. Simi- 
lar states exist at zigzag edges of graphene bilayer^^^^, 
and in other multilayered graphene systems^. Experi- 
mentally, these states were observed in STM experiments 
near monatomic steps on a graphite surface^^^^. 

The presence of large number of localized states is im- 
portant for the predicted edge magnetism in graphene 
nanoribbon^^, a topic that has recently seen renewed 
interest in the context of graphene spintronicJ^^M^. 
Apart from edge magnetism, interacting edge states may 
also result in other correlated ground state j^^tt24| 

Edge states also play a role in confined geometries, 
when the edge to area ratio is large enough so that the 
electronic properties of the boundary may become dom- 
inant. One example for such a geometry are graphene 
quantum dots that have been under intense experimental 
study recentl}*^^*^, with quantum dot sizes in the range 
from a few tens of nanometers to micrometers. Another 
example are antidot arrays that have been subject of sev- 
eral theoretical studies^^^^^ and have also been realized 
exp crime nt all}EM^. 

Different edges have been observed in graphitj ^-*^^ ! ^^ * 
and graphenpSffil. In particular, the existence of bound- 
aries with a long-range crystalline order in exfoliated 
graphene has been questioned^"^. In addition, the exis- 
tence of unsaturated dangling bonds at edges makes them 
reactive, and it is unclear how they are passivated"^^ 
Hence, it is likely that graphene edges are perturbed and 
that the presence of edge distortions has to be taken into 



account. 

The aim of our paper is to show that edge states can be 
expected in realistic disordered quantum dots. We also 
analyse the particular properties of edge states such as 
their number and compressibility. We start the analysis 
in Section [IT] by using the theory of Ref. 9 for a relation 
between the number of edge states per unit length of a 
smooth boundary (see Fig. ^) and the angle the bound- 
ary makes with respect to the crystallographic axis. We 
extend the earlier results by calculating the correction to 
the edge states number coming from the edge roughness 
(Fig. [l]^). Having the total number of edge states and 
their momentum distribution we apply perturbation the- 
ory to see how confinement energy and particle-hole sym- 
metry breaking terms in the Hamiltonian shift the edge 
states from zero energy. Confinement energy spreads a 
delta function-like peak in the density of states into a hy- 
perbolic one. In contrast, particle- hole symmetry break- 
ing terms spread the edge states nearly homogeneously 
over a band of finite width. For realistic dot sizes around 
tens of nanometers we find the latter to be more impor- 
tant. 



In Section HI we perform numerical simulations on 
quantum dots of experimentally relevant sizes. These 
numerical calculations confirm our analytic results. We 
also study the magnetic field dependence of edge states 
in quantum dots. Whereas magnetic field spectroscopy 
of energy levels has up to now mainly been a useful tool 
to probe bulk states in graphene quantum dots^^^^^^, 
we show how to employ this technique also to identify 
edge states. In addition, we study the level statistics of 
edge states. 

Finally in Section [IV| we calculate an upper bound on 
the magnetic moment of a graphene dot due to edge state 
polarization. We also give an upper bound on the relative 
weight of the edge states with respect to the bulk states. 
By knowing the magnitude of additional compressibility 
due to the edge states we estimate parameters of an an- 
tidot lattice in which edge states would be visible in SET 
experiments. 

We conclude in Section [Vl 
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Note that Eq. Q only gives the local density of low 
energy edge states. It should not be confused with Lieb's 
t heor em p2 which connects the number of states with ex- 
actly zero energy with the difference in the number of 
sublattice sites in a bipartitte sublattice. Lieb's theorem 
was applied to graphene in Refs. |3T1 [33l and l52l and for 
a disordered quantum dot geometry it predicts^ number 
of zero-energy modes ~ \flj. Our analysis shows that 
there will be ~ L low energy edge states, although most 
of them do not lie at exactly zero energy. Hence, there is 
no contradiction with Lieb's theorem. 



FIG. 1. A graphene quantum dot. The excess density of states 
due to edge states is shown in a color plot^^, as calculated for 
a quantum dot with a smooth boundary and no particle-hole 
symmetry breaking perturbations (a). In general, edge states 
are present both near a smooth boundary (b) and a boundary 
with short range disorder (c). 



II. ANALYTICAL CALCULATION OF THE 
EDGE STATES DENSITY 

A. Number of edge states 

The density of edge states per unit length was calcu- 
lated for a smooth edge in Ref. 9 
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with — 7r/6 < ^ < 7r/6 the angle boundary makes with a 
nearest armchair direction, and a the lattice constant. 
This expression is valid on the scales larger than the 
boundary roughness scale and another scale depen- 
dent on boundary structure. For most boundary orien- 
tations, except the ones very close to armchair direction 
~ a. Approximating the dot shape by a circle, and 
integrating Eq. [T] along the whole perimeter of the dot, 
we get 
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with L circumference of the dot and a the lattice con- 
stant. This density of states is the difference between 
total density of waves evanescent away from the bound- 
ary and the number of conditions the wave function must 
satisfy on the selected sublattice (see Ref. 9 for a more 
detailed description). If a small fraction a of random 
outermost atoms of the smooth edge oriented at angle 
with armchair direction is etched, the number of con- 
ditions for the wave function on the minority sublattice 
increases by 



SN = 2a sin ( 
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This leads to the reduction of the number of the edge 
states near an edge with atomic scale disorder: 



N' = N{l-2a). 



B. Edge state dispersion 

There are two different mechanisms which give finite 
energy to otherwise zero energy edge states: the overlap 
between edge states on different sublattices, and terms 
breaking sublattice symmetry at the edge. The disper- 
sion resulting from these perturbations can be calculated 
by applying degenerate perturbation theory, acting on 
the wave functions ^/^n, belonging exclusively to A or 
B sublattice. The long wavelength part of these wave 
functions is defined by the conformal invariance of Dirac 
equation, so they can be approximated as plane waves 
belonging to one of the six facets of the dot with well- 
defined boundary condition, extended along the facet and 
decaying into the bulk. These wave functions have lon- 
gitudinal momenta 
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approximately equally spaced due to phase space argu- 
ments. 

We first estimate the energy dispersion due to edge 
state overlap, or in other words by finite size effects. 
Particle-hole symmetry prevents coupling between states 
on the same sublattice, so the dispersion of edge states in 
a finite system can be calculated from the matrix element 
between the edge states on different sublattices. These 
states are separated from each other by a distance of 
an order of the dot radius R and their decay length away 
from the boundary is proportional to difference k between 
their momentum and the momentum of the nearest Dirac 
point (Dirac momentum), so the energy is 
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where vy is the Fermi velocity and we set = 1. We note 
that Eq. [6] is very similar to the energy of edge states 
in zigzag nanoribbons^^. Substituting the value of mo- 
mentum of the edge states from Eq. ([5| into Eq. ([6| we 
calculate the density of edge states per unit energy 
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The atoms passivating the edge perturb the 7r-orbitals 
of carbon atoms to which they are bound. This in- 
teraction breaks the effective electron-hole symmetry of 
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graphene around the Dirac point. Next-nearest neigh- 
bor hopp ing is breaking this symmetry at the edges as 
welP^^, and it was shown to be equivalent to the edge 
potential. The dispersion of the edge states near a 
zigzag edge due to these two perturbations is 

£; = (Ae-tO[cos(K)-l/2], 27r/3 < K < 47r/3, (8) 

where K is the full momentum of the edge state, Ae 
is the average strength of the edge potential and t' the 
next-nearest neighbor hopping strength. Despite it is 
not straightforward to generalize this equation to an ar- 
bitrary orientation of the edge, the general effect of the 
electron-hole symmetry breaking terms is to smear the 
zero energy peak in the density of states into a band 
between energies of approximately and Eq = /S.e — t' 
for the most localized states, while the more extended 
states are near the Dirac energy. The one dimensional 
van Hove singularity in the density of states E = Eq 
will be smeared out, due to the presence of a minimum 
decay length of the edge states when the orientation of 
the boundary is not exactly zigzag.^^ 

The energy due to finite size of the dot given by Eq. (|6| 
is at best of an order of E ^ / R ^ ta/ R. It is less than 
tens of millivolts for dots above lOnm size. On the other 
hand the energy due to the edge potentials and next- 
nearest neighbor hopping (Eq. ([8|) is likely to be around 
hundreds of millivolts. Accordingly in realistic dots with 
edge potentials and next-nearest neighbor hopping term 
edge states occupy the band between the Dirac point and 
£^0 with approximately constant density 

R 

Pedge = c(l-2a)^^, (9) 
with c = 8-4^3^ 1. 



III. NUMERICAL RESULTS 

In order to confirm the analytical results of the previ- 
ous sections we have performed numerical simulations of 
the energy spectrum of graphene quantum dots with sizes 
relevant for experiments. In the following we present re- 
sults for a quantum dot with the shape of a deformed 
circle'^ (c.f. Figjl]), characterized by an average radius 
R. Although we focus on a particular quantum dot here, 
we have found through numerical studies that the char- 
acteristic features of our results are independent from the 
details of the dot shape. 

The numerical simulations are based on a tight-binding 
model of graphene with Hamiltonian 

H = -^tijclcj + h.c. (10) 

where the hopping tij = t for nearest neighbors and tij = 
t' for next-nearest neighbors^. The effects of a magnetic 



field are incorporated through the Peierls phase a^ 

tij tij X exp J dsA(x)^ , (11) 

where and Xj are the positions of atom i and j, re- 
spectively, and A(x) is the magnetic vector potential. 

The quantum dots are constructed by "cutting" the 
desired shape out of the hexagonal graphene grid. For 
a shape that is smooth on the length scale of the lat- 
tice constant as considered here, this results in edges 
with a locally well-defined orientation {smooth edges^ see 
Fig. [ijb)). In order to account for edge disorder on the 
lattice scale {rough edges^ see Fig. [ijc)), we adopt the 
disorder model introduced in Ref. 59: Starting from the 
smooth edge, atoms at the boundary are removed ran- 
domly with probability p, with dangling bonds removed 
after each pass. This procedure is repeated A^sweep times. 

The energy spectrum of the dot tight-binding Hamil- 
tonian is calculated numerically using standard direct 
eigenvalue algorithm^ and matrix bandwidth reduction 
techniques'^ if a large part of the spectrum is needed. 
In contrast, if only a few eigenvalues and -vectors are 
sought, we apply an iterative technique'^ in shift-and- 
invert mode^. 



A. Systems with electron-hole symmetry 

We first focus on the electron- hole symmetric case, 
i.e. t' = {) and the absence of potentials. Fig. [2|a) shows 
the number of states N{E) per energy interval /S.E for 
dots with smooth and rough edges. We can clearly iden- 
tify the edge states close to £^ = and the linearly in- 
creasing bulk density of states. Approximating the cir- 
cumference of the dot as L ^ 27ri?, Eq. ([2| predicts 
^ 170 edge states for a quantum dot with a smooth 
edge, which is in very good agreement with A^ = 169 ± 6 
edge states obtained from the numerical simulation by 
summing over the three central bins, where the number 
of states differs noticeably from the linear bulk density 
of states. The number of edge states N' in the dot with 
atomic scale disorder can be estimated from Eq. (|3| by 
approximating a ^ pA^sweep yielding A"' ^ 0.5A^ for the 
disorder parameters used in the simulation (A^sweep = 5, 
p = 0.05), again in good agreement with the numerical 
simulations. 

In order to examine the behavior of the edge state den- 
sity of states in more detailed close to £^ = 0, we estimate 
the density of states numerically as 

p{{Ei+, + Ei)/2) = ^ , (12) 

where Ei is the energy of the i-th state in the dot. 
Fig.|2|b) shows the numerically computed p{E) for quan- 
tum dots with smooth and rough edges. As predicted in 
Eq. ([7|), we find a l/£^-dependence close to E = 0; quite 
remarkably, we find an excellent agreement with this scal- 
ing for more than ten orders of magnitude. The clustering 
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FIG. 2. Electronic states in a graphene quantum dot close to 
the Dirac point. The graphene quantum dot has the shape 
of a deformed circle (see Fig. [l] and footnote 57) with R = 
160a ~ 40 nm, and we consider both smooth and rough edges 
as shown in Fig. [ijb) and (c) respectively. The parameters 



for the edge disorder are A^s^ 



5 and p = 0.05 (see the 



main text for a discussion of the edge disorder model), (a) 
Number of states per energy interval AE for a quantum dot 
with smooth (black lines) and rough edges (red lines), with 
= 0. 4t/6 1. (b) Density of states estimated numerically 
from Eq. ( 12 ) for a quantum dot with smooth (black symbols) 
and rough edges (red symbols). For comparison, the blue 
dashed line shows a l/£'-dependence. 



of data points at p{E) = 10^^ is due to the finite pre- 
cision in the numerical calculations. It should be noted 
that we found this remarkable agreement with theoretical 
predictions without averaging over an energy window or 
different dot shapes, implying that the spectrum of edge 
states is highly non-random even in a quantum dot with 



random shape. We come back to this point in Sec. HID 



B. Broken electron-hole symmetry 

Next we focus on perturbations breaking the electron- 
hole symmetry. For this we consider a finite next-nearest 
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FIG. 3. Number of states (black lines) per energy interval 
for a graphene quantum dot with smooth (left panels) 
and rough (right panels) edges. We show results for situations 
with broken electron-hole symmetry: (a) finite next-nearest 
neighbor hopping and no edge potential {t' = O.lt and Uo = 0) 
and (b,c) finite next-nearest neighbor hopping including an 
edge potential (t' = O.lt, with (b) Pedge = 0.25 and Uo = 0.2t, 
and (c) Pedge = 1 and Uo = O.lt). The remaining parameters 
are as in Fig. [2] The blue dashed lines show the number of 
bulk states A^buik estimated from the linear density of states 
of the Dirac dispersion Eq. (13). 



neighbor hopping t' uniformly within the quantum dot, 
as well as a random potential at the quantum dot edge, 
where an energy Uq is assigned to edge atoms with prob- 
ability Pedge- 

Fig. (sja) shows the number of states per energy win- 
dow AE for finite t', but in the absence of an edge po- 
tential. In order to identify the edge states properly, we 
compare the numerical data including the edge states to 
the number of bulk states estimated from the linear Dirac 
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density of states 
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approximating the area of the quantum dot as A = iiR^ . 
The excess edge state density of states can be clearly 
identified, both in the case of smooth and rough edges. 
The bulk density of states close to £^ = is unaffected 
by a finite t\ the effect of electron-hole asymmetry on 
the bulk states only shows for energies \E\ > O.lt. The 
central edge state peak observed for f = (c.f. Figj2| 
is broadened and shifted towards the hole side, but the 
total number of edge states remains unchanged from the 

= case. The excess density due to the edge states 
is approximately constant in the energy range between 

= —O.lt and 0, in accordance with the prediction from 
Eq. (|9|. As before, atomic scale edge disorder only results 
in a reduction of the total number of edge states. 

The presence of an additional edge potential changes 
the energy range of the edge states. In Fig.jsj^b) we show 
results for an average edge potential Ae = O.Obt. Cor- 
respondingly, the majority of the edge states occupies 
uniformly an energy window between Ae — t' = — 0.05t 
and 0. A few states can still be found beyond this en- 
ergy window, as the randomness of the edge potential 



has been neglected in the arguments of Section |IIB| In- 
stead, if the edge potential is uniform, the dispersion of 
the edge state due to next-nearest neighbor hopping can 
be cancelled exactly by Ae = — t', as shown in Fig. [sj^c). 
This particular example strikingly shows the equivalence 
of next-nearest neighbor hopping and an edge potential, 
as predicted in Ref. 56, 

The narrowing of the energy band width occupied by 
the edge state due to an edge potential may also be a 
possible explanation (amongst others^^) for the fact that 
STM measurements on zigzag graphene edges found a 
peak in the density o f stat es only a few tens of meV 
below the Dirac point far less than expected from 
estimated values of the next-nearest neighbor hopping. 



Broken time- reversal symmetry: 
magnetic field 



Finite 



We now consider the effects of a finite magnetic field on 
the edge state energies. The evolution of edge states in a 
magnetic field has been studied theoretically for s pecial 
geometries and a particle-hole symmetric spectrunP^'^. 
Recently, the magnetic field dependence of the energy 
levels in a graphene quantum dot has been also been sub- 
ject to an experimental investigation^^. However, in the 
theoretical calculations used to interprete these experi- 
ments the graphene edge states were excluded artificially. 
As we show below, the presence of edge states results in 
a much richer magnetic field dependence of energy lev- 
els in a graphene dot, in particular when particle-hole 
symmetry is broken. 
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FIG. 4. Magnetic field dependence of the energy levels (black 
lines) in a desymmetrized quantum dot with R = 100a (de- 
formed circle as shown in Fig. [l] c.f. footnote 57). The par- 
ticipation ratio p of the states is color-encoded, with the most 
strongly localized states in red. The blue dashed lines indi- 
cate the energy of the n — 0, ±1 Landau levels of graphene. 
The calculations includes finite next-nearest neighbor hopping 
t' — O.lt and a random edge potential with Pedge = 0.25 and 
Uq = 0.2t. 



In Fig. |4]we show the numerically calculated magnetic 
field dependence of the energy levels in a graphene quan- 
tum dot close to the Dirac point, for finite t' and edge 
potential. In order to distinguish between edge and bulk 
states, we also plot the participation ratic 



P 



(14) 



where the index i runs over atomic sites and N denotes 
the total number of atoms in the dot. The participation 
ratio p can be interpreted as the fraction of atoms occu- 
pied by an electron for a given energy level. Thus, p ^ 1 
for extended states {p ~ 0.3 — 0.4 in quantum dots) and 
p <C 1 for localized edge states {p ^ 10~^ — 10~^). 

Instead of a uniform fiow of energy levels towards the 
n = Landau level as calculated in Ref. 28 , we observe that 
the most strongly localized states only show a very weak 
magnetic field dependence (apart from avoided cross- 
ings), leading to a far richer energy spectrum. Note that 
this effect is most prominent on the hole side of the spec- 
trum where the majority of the edge states reside, as can 
be simply seen by comparing the number of states for 
^ > and E <{). This weak magnetic field dependence 
of the localized edge states can be understood from the 
fact that bulk states start to be affected by the magnetic 
field when the cyclotron radius becomes comparable to 
the dot size, whereas edge state energies are expected to 
only change significantly when the cyclotron radius be- 
comes comparable to the edge state decay length which 
is much smaller than the dot dimensions. 

Note that this type of behavior is similar to the mag- 
netic field dependence of the low-energy spectrum of 
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graphene in the presence of lattice vacancies!^ In fact, 
such vacancies can be considered as internal edges and 
also carry a localized state. 

Hence, magnetic field independent energy levels are 
characteristic for localized (edge) states. In the light of 
this observation, it would be very interesting to see if 
experiments can identify such states, which would be a 
strong indication for the presence of such states. 



D. Level statistics of edge states 

The bulk states of chaotic graphene quantum dots con- 
fined by lattice termination have been shown to follow 
the level statistics of the Gaussian orthogonal ensem- 
ble (GO E), as expected for a system with time-reversal 
symmetr}^^!^ (scattering at the quantum dot boundary 
mixes the K and i^'-valley). The edge states however 
are tied to the boundary of the quantum dot only, and 
should not necessarily follow the same level statistics as 
the extended states. Instead, being localized states they 
are rather expected to follow Poisson statistics, as has 
also been noted in Ref. |70l but not been demonstrated 
explicitly. 

To check these expectations we have studied the level 
spacing distribution of edge states in quantum dots. For 
this purpose we have identified edge states using the par- 
ticipation ratio and worked with the edge state spec- 
trum alone. This spectrum has been unfolded''^ using 
the average density of states and scaled to an average 
level spacing of unity. The distribution P{S) of the 
nearest-neighbor level spacings S in the unfolded spec- 
trum is then normalized such that / P{S)dS = 1 and 
j SP{S)dS=l. 

Fig. [5] shows the level spacing distributions for the 
electron-hole symmetric case (t' = 0) and for broken 
electron-hole symmetry {t^ = O.lt). Surprisingly, the 
edge states follow the GOE statistics if f = 0. Only 
if a finite is included, they exhibit a statistics close to 
Poisson. These classifications are additionally corrobo- 
rated by the integrated level spacing distributions shown 
in the inset of Fig. |5j 

This striking difference in level statistics can be ex- 
plained by the different nature of the wave functions: 
The graphene Hamiltonian exhibits a chiral symmetry for 
t' = that results in an equal occupation probability of 
sublattice A and B for every individual wave function.!^ 
Since the edge wave function at a certain type of zigzag 
edge is nonzero only on one sublattice, every eigenstate 
for t' = must also occupy another part of the boundary 
of the opposite kind, as illustrated in Fig. [6) This leads to 
an artificial long-range coupling between edge states and 
thus to level repulsion, resulting finally in GOE statis- 
tics. If this chiral symmetry is broken, for example by 
next-nearest neighbor hopping,'^ edge state wave func- 
tions may be localized at a single edge only (Fig. [6]^ a)). 
Whereas edge states localized at the same part of the 
boundary still may feel level repulsion, parts that are 




FIG. 5. Level spacing distributions for quantum dots with 
smooth edges for = (solid red curve) and t' = O.lt 
(solid black curve), together with the theoretical predictions 
for Poisson statistics (dashed line), the Gaussian orthogonal 
ensemble (dash-dotted line), and the Gaussian unitary en- 
semble (dotted line). The inset shows details the integrated 
level spacing distribution for small level spacings S (same 
line colors and -types as the main plot). The level distri- 
bution statistics has been obtained by averaging individual 
level distributions from 100 quantum dots similar to the type 
given in footnote i57i with average radius R = 160a. A state 
has been identified as an edge state, if its participation ratio 
Pi < OmJ^ For t' = we have also omitted all states with 
an energy smaller than the numerical precision. 



t' = 



t' = o.lt 
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FIG. 6. Color plotP of wave function density in a graphene 
quantum dot (shape as described in footnote I57|) for the e- 
h symmetric case (t^ = 0, left column) and for broken e-h 
symmetry (t' — O.lt, right column) on the examples of a 
mode that is (a) strongly decaying and (b) slowly decaying 
into the bulk. Note that for presentation purposes we have 
chosen a rather small dot (R = 30a), but the behavior does 
not change qualitatively for larger dots. 
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further away may only interact via hybridization with 
bulk states which typically happens for edges states de- 
caying further into the bulk, as seen in Fig.loFb). For the 
type of quantum dots under consideration (Fig. [T]), this 
results in six approximately independent series of energy 
levels, and hence an approximate Poisson statistics. 

A finite next-nearest neighbor hopping t' (or another 
chiral symmetry breaking term) thus does not only 
change properties of the edge states quantitatively, but 
leads to a striking, qualitatively different level statistics. 



point. From Eqs. Q and (13), the average ratio between 



edge and bulk states in this energy range is 

c{l-2a) 



Nhulk 



ElaR 



(17) 



This gives, for a diameter of lOOnm and Eq = 0.3eV an 
upper bound of Nedge/Nbuik ^1/2. 



Detection in antidot lattices 



IV. DISCUSSION AND PHYSICAL 
IMPLICATIONS 

A. Formation of magnetic moments at the edges 

An extensively discussed topic in the graphene lit- 
erature is the formation of localized moments at 
boundarie4iIii^l™t^lZ6llll. 

The previous analysis al- 
lows us to set approximate bounds on the maximum mag- 
netic moment in a graphene quantum dot. 

The interaction energy between two electrons of oppo- 
site spin in a boundary state of area k^^ x R ^ a x R 



E., 



e\ (R 
^ log - 
R \ a 



(15) 



where e is the electronic charge. States with energies 
\ei — Ef\ ^ will be spin polarized. Since the density 
of edge states is nearly constant and given by Eq. (|9]), 
the position of the Fermi level is not relevant. Using 
the density of states given in Eq. ([9|, we obtain for the 
number of spins in a quantum dot: 



R 

spins ~ EeePedge = c(l - 2a) ^-^ log ' 



20(1 - 2a) log 



R 



(16) 



where for last estimate we took Eq = 0.3 eV. The maxi- 
mal number of polarized spins depends only logarithmi- 
cally on the size of the dot. 

In general, the states at the edge of a quantum dot 
will belong to one of the two sublattices with equal prob- 
ability. States localized at different sublattices interact 
antiferromagneticall}'^^ . If we neglect this interaction, 
we expect a maximum magnetic moment comparable 
with A^spins- When the antiferromagnetic interaction con- 
tributes to the formation of the total magnetic moment, 
its value will be proportional to the number of uncom- 
pensated sites at the edges, which will scale as y/N^^i^. 



B. Fraction of edge states 

Our results suggests that edge and bulk states can co- 
exist in a range of energy of order Eq near the Dirac 



A conclusive way of detecting the existence of edge 
states can be the measurement of their contribution to 
the electronic compressibility. It is hard to detect the 
edge states in a single quantum dot because the ground 
state properties are dominated by the charging energy. 
Also, the contribution of edge states to the density of 
states in most large-scale samples will be negligible com- 
pared to the bulk contribution. However it is possible 
to circumvent both problems in antidot lattices. The 
Coulomb energy does not play a role in this case due to 
absence of confinement. On the other hand, the existence 
of multiple antidots allows us to reach a large edge-area 
ratio. To estimate whether it is possible to detect edge 
states, we use the value of minimal compressibility (or 
the minimal density of states) of bulk graphene Ref. (79] 



dn 



3 X IQ-^^meVcm^ 



(18) 



and we assume that the width of the band of edge states 
is around Eq ^ 0.3 eV. 

We consider an antidot lattice with antidot size L of 
the same order of magnitude as the antidot spacing. Us- 
ing the analysis in the previous section, the density of 
states per unit area associated to the edge states is: 



N-UE)^EoaL 



(19) 



Comparing this expression with eq. [18) and using Eq ^ 
0.3 eV, we find that the contribution from the edge 
states is comparable to the bulk inverse compressibility 
for L < Ijum. Hence, the additional density of states 
near the edge will be visible in compressibility measure- 
ments using a single electron transistor (SET)since the 
size of the SET tip is around 100 nanometers.^ Our re- 
sults may be the reason of p-doping observed in antidot 
lattices experiment ally. ^-^^ 



V. CONCLUSIONS 

We have analyzed generic properties of the electronic 
spectrum in graphene quantum dots. We find that some 
of the electronic states will be localized at the edges and 
form a narrow band. The density of states in this band 
is (X 1/E in graphene dot without electron- hole symme- 
try breaking perturbations. In presence of such pertur- 
bations, the density of the edge states is approximately 
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constant and scales as R/aEo, where R is the dot radius, 
a is the lattice constant, and Eq is an energy scale which 
describes the edge potentials and next-nearest neighbor 
hopping. 

If chiral symmetry is present, the edge states expe- 
rience strong level repulsion and are described by the 
Gaussian orthogonal ensemble. Chiral symmetry break- 
ing terms (such as next-nearest neighbor hopping) how- 
ever lift this spurious level repulsion leading to the Pois- 
sonian statistics expected for localized states. In con- 
trast, extended states will be described by the orthogonal 
or unitary ensembles, depending on t he str ength of the 
intervalley scattering at the boundarie^^^. 

Having an analytical model for the edge states allows 
us to estimate the maximum spin polarization due to the 
presence of edge states. We predict that the additional 
density of states due to edge states will be visible in SET 



experiments. Effect of edge states on transport in quan- 
tum dots and more detailed investigation of interaction 
effects remains a direction for further research. 
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